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Abstract The supernova yields of r-process elements 
are obtained as a function of the mass of their progenitor 
stars from the abundance patterns of extremely metal- 
poor stars on the left-side [Ba/Mg]-[Mg/H] boundary 
with a procedure proposed by Tsujimoto and Shigeyama. 
The ejected masses of r-process elements associated with 
stars of progenitor mass M ms < 18M Q are infertile sources 
and the SNe II with 20M© < M ms < 40M©are the dom- 
inant source of r-process nucleosynthesis in the Galaxy. 
The ratio of these stars 2OM < M ms < 40M© with 
compared to the all massive stars is about ^18%. In this 
paper, we present a simple model that describes a star's 
[r/Fe] in terms of the nucleosynthesis yields of r-process 
elements and the number of SN II explosions. Combined 
the r-process yields obtained by our procedure with the 
scatter model of the Galactic halo, the observed abun- 
dance patterns of the metal-poor stars can be well repro- 
duced. 
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1 Introduction 

Elemental abundances in metal-poor Galactic halo stars 
are providing evidence of the earliest Galactic nuclc- 
osythesis history and clues about the identities of the 
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first stellar generations, the progenitors (or predecessors) 
of the halo stars. The sample of identified field stars 
with [Fe/H]<-2.5 has increased by more than an order 
of magnitude in the last decade. Because the elements 
in the atmospheres of these stars have been produced 
in a small number of nucleosynthetic events, abundance 
determinations can provide direct tests of model yields 
from different nuclear processes. 

The elements heavier than the iron peak are made 
through neutron capture via two principal processes: the 
r-process and the s-process (Burbidge et al., 1957). The 
r-process (for rapid process) occurs when neutrons are 
added much more rapidly than the (3 decay times of the 
relevant nuclei. The site or sites of the r-process are not 
known, although suggestions include the ^-driven wind 
of Type II SNe (e.g., Woosly and Hoffman, 1992; Woosly 
et al., 1994) and the mergers of neutron stars (e.g., Lat- 
timer and Schramm, 1974; Rosswog et al., 2000). 

Particular attention to the Galactic evolution of ele- 
ments produced by neutron-capture nucleosynthesis was 
given by Mathews et al., (1992), Pagel and Tautvaisiene 
(1997), and more recently by Travaglio et al. (1999). 
These authors adopted the standard approach to Galac- 
tic chemical evolution, assuming that stars form from a 
chemically homogeneous medium at a continuous rate. 
A more realistic model for the chemistry and dynamics 
of the gas is needed in order to investigate the earliest 
phases of halo evolution. 

Recently, several chemical evolution models statis- 
tics on the build-up of chemical elements in the early 
Galaxy (e.g., Tsujimoto et al., 1999; Argast et al., 2000; 
Oey, 2000; Travaglio et al., 2001; Fields et al., 2002); 
Tsujimoto et al., (1999) provided an explanation for the 
spread of Eu observed in the oldest halo stars in the con- 
text of a model of supernova-induced star formation. As- 
suming that r-process nucleosynthesis sites would most 
likely be identified with Type II supernovae (SNe II), 
Mathews et al. (1992) suggested a mass range of M ms =7— 
8Mq for the site, while Travaglio et al. (1999) supported 
a somewhat higher mass range, M ms =8— 10M Q . Cescutti 
et al. (2005) concluded that the Eu should originate as 
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an r-process element in stars with masses in the range 
10-30M Q . 

Recent observations and the analyses imply that the 
abundance pattern of an extremely metal-deficient star 
with [Fe/H]<-2.5 may retain information of a preced- 
ing single supernova (SN) event or at most a few SNe 
(Mcwilliam et al., 1995; Ryan et al., 1996). Tsujimoto 
and Shigeyama (1998) have also shown that the mass of 
r-process elements ejected by each SN II as a function 
of progenitor mass at the main sequence (M ms ) can be 
derived from the observed [Ba/Mg]-[Mg/H] trend com- 
bined with the [Mg/H]-M ms relation in theoretical SN 
models. But some of extremely metal-deficient stars with 
[Fe/H]<-2.5 could be formed out of gas enriched by sev- 
eral SNe (Fields et al., 2002). 

In this paper, we assume that the stars on the left-side 
[Ba/Mg]-[Mg/H] boundary are made from individual su- 
pernova events. In Section 2 we discuss the r-process 
elements yields based on the model proposed by Tsu- 
jimoto and Shigeyama (1998). In Section 3 the r-process 
elements yields are used to explore the early stages of 
inhomogencous chemical evolution in the Galaxy. Con- 
clusions are given in Section 4. 



2 Production site and yield for r-process 
elements 

Assuming that the stars on the left-side [Ba/Mg]-[Mg/H] 
boundary retain the abundance pattern of a single su- 
pernova and using the similar procedure presented by 
Tsujimoto and Shigeyama (1998), we can calculate the 
r-process elements yields as a function of the initial stel- 
lar mass from theoretical SN models (e.g., Woosley and 
Weaver, 1995; Tsujimoto et al., 1995; Nomoto et al., 
1997). 

In Fig. la, which uses the data of Mcwilliam et al. 
(1995), Mcwilliam (1998), Lai et al. (2004), Honda et al. 
(2004) and Barklem et al. (2005), the [Ba/Mg] values for 
a sample of metal-poor stars are plotted against [Mg/H]. 
We infer yields for Ba element from the left-side bound- 
ary of observed abundances in stars with -4< [Fc/H] <-2.5 
by calculating the chi-squared fit to the data. The fitted 
line is shown in Fig. la, which represents the mass range 
of the r-process site. This gives a relation between the 
metallicity [Mg/H] of stars and the mass Mms of SN II 
progenitor as shown in Fig. lb. Since the observed abun- 
dances in CS22892-052 and other stars strongly enriched 
r-elements ([Eu/Fe]>1.0) may not reflect the composi- 
tion of the ISM from which it formed (Qian and Wasser- 
burg, 2003, Barbuy et al. 2005), we exclude these stars 
in our calculation. 

Using the theoretical nucleosynthesis mass of Mg for 
each Mms (Tsujimoto et al., 1995), the ejected mass 
of element as a function of M ms is derived. The yield 
of Ba obtained is shown in Fig. 2. In addition, Eu is 
a tracer of the r-process, so we investigate the enrich- 
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Fig. 1 (a) The correlation of [Ba/Mg] with [Mg/H] for 
metal-stars; the thick line is the left-side boundary; the open 
diamonds represent data of stars with [Eu/Fe]>1.0. (b)SN II 
progenitor mass M ms is plotted as a function of [Mg/H] in- 
side the shell swept up by the SNR. Dashed lines in the lower 
and upper panels illustrate how we assign M ma to each star 
along the boundary. 



ment of Eu as a representative of r-process elements. 
Because the abundance pattern of heavy neutron cap- 
ture elements (Z>56) for each star is quite similar to 
that of the r-process component in solar-system mate- 
rial (Senden et al., 1994, Cowan et al., 1995), the Eu 
yield can be derived from M e j t Eu — M e j^B a (^; )©and 
also shown in Fig. 2. According to our present models, 
SNe II with M ms = 12 — 20Af Q do not produce signifi- 
cant amounts of r-process elements. It is indicated that 
SNe with M ms > 18M Q are r-process main sites. The de- 
rived mass of Ba synthesized in SNe II is 3 x 1O~ 6 M for 
M ms = 2OM , which is in good agreement with the val- 
ues of Tsujimoto and Shigeyama (2001) and Pastorello 
(2005). We remind the reader that our model is based 
on the observed abundances of the metal-poor stars, so 
the uncertainties of those observations will be involved 
in the model calculations. We note from Fig. la and lb 
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Fig. 2 Ejected masses of Ba and Eu from SNe II as a function 
of the mass of their progenitor. 



that for the massive stars (M> 25M ) the boundary is 
not explicit, so the uncertainties of the r-process yields 
for these stars are larger than those of lower mass stars. 



3 The r-process scatter and the Galactic halo 
evolution 

The very metal-poor stars, presently found in the halo of 
the Galaxy, are believed to have formed at the ear liest 
times, shortly after it became possible for the Universe 
to make stars with sufficiently long main- sequence life- 
times to survive for ~14 Gyr (Senden et al., 2003). The 
chemical compositions of these stars are thus expected 
to reflect a quite small number of nucleosynthesis pro- 
cesses, possibly as small as one (e.g., the stars on the 
left-side boundary in Fig. fa), while the compositions 
of metal-rich stars reflect the cumulative results of the 
various processes that have been in operation during the 
entire history of Galactic chemical evolution. Models for 
r-process nucleosynthesis and for the Galactic chemical 
evolution must account for the scatter of r/Fe between 
halo stars, as well as for the smaller scatter between disk 
stars. We assume, as Fields et al. (2002) did, that the 
observed Pop II r/Fe abundance in each star reflect r- 
process contributions of a few SNe II. The basic idea is 
that the r/Fe scatters arise from mixing among differ- 
ent mass SNe II, which can produce different yields of 
r-process elements. 

To model the scatter of r/Fe, we create a set of halo 
stars and deduce the history-the nuceosynthctic ancestry 
-of each. We assume that a single supernova event mixes 
its metals and r-elements over a region of mass M and 
analyze the mixing history of the interstellar medium 
(ISM). For the first generation of SN, the r-elements and 



iron yield of in the ith region is m r i(i) and m Fe i(i) re- 
spectively. Since a given parcel of ISM gas can be en- 
riched to different degrees by the ancestors it had, the 
new region (with mass M also) polluted by the successive 
generation of SN comprises material originating from J 
old regions. The accumulated r-elements and iron yields 
of ith new region after the Nth generation of SN are 



m rN (i) = m' r 



(JV-l),i 



+ m rN 



m FeN (i) = m'p^x^i + m FeN 



(1) 
(2) 



where the reduced yields of the former (N-l) gen- 
erations SNe contributed to the new region arc 



m 



r(JV-l), 



j = l- j,(iV-l)"V(JV-l)(fc.j) andm fe(JV-l),i = 

i- i ~'j,(N-i) m Fe(N-i)(kj). Here we randomly take J 
group yields from the (N-l)th generation regions, which 
are tn r (N-i){kj) and vn Fe i^_i^{kj) (j = l> • ■ • >J) respec- 
tively. Random parameters Pj^N-i) (j=lj • • • >J) are t ne 
mass fractions (statistical weight) came from J old re- 
gions to form the zth new region which will be polluted 
by Nth generation SN. For P,-,(jv-i), Ej=i p j = l - 

To simplify notation, we will define the scale r/Fe 
ratio to be 



R 



-/Fe 



(r/Fe) e 



(3) 



which implies that [r/Fe]=logR When taking the yield 
of SN II with M ms = 20M© and m Fe = 0.07M Q , we 
can obtain a value of R«50, which is consistent with the 
value of Fields et al. (2002). We assume that each of the 
halo stars we create incorporates gas cloud which has 
been enriched by some number N of supernova: this is 
the number of supernova ancestors for the star. Let the 
total number of supernova events be N, we have that 



m rN (i)/A r 
R N,i = 7T71 — /( r / Fe )& 



mFeN(i)/A 



Fe 



(4) 



At the very earliest time, when single supernova events 
really contaminate a given star, the [r/Fe] scatter only 
populates the extremes; at later time, mixing becomes 
efficient gradually, and then the scatter decreases. So 
the [r/Fe] becomes a tracer of the inhomogeneity of the 
halo. Several groups (Argast et al., 2000; Tsujimoto et 
al., f999; Oey, 2000) have used similar arguments 
to motivate detailed models that explain the observed 
[r/Fe] (specifically, [Eu/Fe]) scatter in terms of r-process 
nucleosynthesis and an inhomogencous chemical evolu- 
tion of the Galactic halo. In our model, the Eu yield is 
given in Fig. 2 and iron yield is adopted as follows: 
Case A: M Fe = O.O7M 

Case B: M Fe is a function of progenitor mass (Tsujimoto 
et al., 1995). 

The metallicity determines the total number of su- 
pernova ancestors via 



[Fe/H] N>i = logX Fe * - logX FeQ = log 



m FeN (i) 
M x X FeQ 



(5) 
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Fig. 3 Correlation of [Eu/Fe] with [Fe/H]. The symbols rep- 
resent the data taken from (open circles), Gilory et al. (1988) 
(+ center squares) , Hartmann and Gehen (1988) (open 
hexagons) ,Zhao and Magain (1990), Magain (1989) (open 
squares) , Da Siler et al.(1990) (open down triangles) , Pe- 
terson et al. (1990) (pluses), Edvardsson et al. (1993) (cross 
center squares) ,Gratton and Sneden (1994) (open up tri- 
angles), Woolf et al. (1995) (open stars) , Mcwilliam et al. 
(1995), Mcwilliam (1998) (open circles), Shetrone (1996) (+ 
center diamonds), Jehin et al. (1999) (open diamonds), Bur- 
ns et al. (2000) (- center squares), Cowan et al. (2002) (star), 
Honda et al. (2004) (open pentagons), Barklem et al. (2005) 
(crosses). The dots represent the Case A result of our model 

where Xp e Q and Xp e are the mass fraction of Fe in solar 
system and the polluted region respectively. In the early 
Galaxy, the gas composed of only hydrogen and helium 
with their ratio X# :A^ff e =0.75:0.25, so 

M = M sw /0.75 (6) 

where M$w is the mass of hydrogen swept up by an su- 
pernova remnant (Shigeyama and Tsujimoto 1998). On 
one hand, when J is large enough, equations (4) and (5) 
will give the form of well-mixed chemical evolution of 
the Galactic halo which could not explain the scatter of 
[Eu/Fe] observed for metal-poor stars, on the other hand, 
when J=l, these equations will give the result that all 
the Fe and r-elements ejected from N supernovae may 
pollute only one supernova remnant and this is not what 
actually occurs. In fact, the J must be very small, so we 
take 3=2 in this work. 

It is of interest to consider making the mass of an- 
cestor a random variable. The initial mass distribution 
of supernova progenitor is generated according to the 
formula of Eggleton et al. (1989), 



(l-Xy>- 75 + 0.032(1 -X) - 25 



Fig. 4 Correlation of [Eu/Fe] with [Fe/H]. The symbols rep- 
resent the data taken from the same refs of Fig. 3. The dots 
represent the Case B result of our model. 

where X is a random number uniformly distributed be- 
tween and 1, which leads to a mass function similar 
to that of Miller and Scalo (1979). Results for a Monte 
Carlo simulation of a stellar population appear in Fig. 
3 and Fig. 4 for Case A and Case B respectively. The 
degree of scatter increases with de creasing metallicity 
because of counting statistics, so that the lowest metal- 
licity events record the nucleosyn thesis of a few events. 
Using the abundance pattern of heavy neutron-capture 
elements (Z>56) is quite similar to that of the r-process 
component in solar-system material (Sneden et al., 1994; 
Cowan et al., 1995), Figs. 5 and 6 show the scatters of 
[Ba/Fe], [Ce/Fe], [La/Fe], [Nd/Fe], [Pr/Fe], [Sm/Fe] for 
Case A and B respectively. A comparison with the ob- 
served data in Figs. 3-6 shows that the model gives a 
good fit to the available data. We thus conclude that 
the observed scatter in [r/Fe] can be understood by our 
model. 



4 Conclusions 

Assuming that the stars on the left-side [Ba/Mg]-[Mg/H] 
boundary retain the abundance pattern of a single super- 
nova, we have calculated the r-process elements yields 
as a function of the initial stellar mass from theoret- 
ical SN models. The scatter of [r/Fe] in halo strongly 
suggests that SNe II associated with stars of progenitor 
mass M ms < 18M Q are infertile sources for the produc- 
tion of r-process elements, which is in agreement with the 
result of Tsujimoto et al. (2000). We conclude that SNe 
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Fig. 5 Correlation of [X/Fe] with [Fe/H]. The symbols represent the data taken from the same refs of Fig. 3. The dots 
represent the Case A result of our model. Where X represents successionally Ba, Ce, La, Nd, Pr, Sm. 



II with 20M Q < M ms < 40Mq are the dominant source 
of r-process nucleosynthesis in the Galaxy. The effect on 
the SNe II which produces the high r-process yield could 
be significant. Assuming a Salpeter IMF with a=-2.35, 
the ratio of these stars compared to the all massive stars 
is about ~18%, which is higher than the value of case A 
in Fields, Truran, and Cowan (2002). In fact, the mean 
value of [r/Fc] obtained in our work is also higher than 
the value of R in Fields et al. (2002). 

We have presented an approach to understand the 
scatter in heavy r-process-to-iron ratio in metal-poor halo 
stars and used a stochastic description of progenitor mass 
SNe II that contribution to the heavy r-process and iron 
abundances in each halo star. The random star-to-star 
variations in nucleosynthctic ancestry lead to scatter in 
[r/Fe]. The models we present are all successful in repro- 
ducing the scatter in the available data, which go down 
to about [Fe/H]«-3. In conclusion, the ratios of [r/Fe] 



with the metallicity are twofold. One is the abundance 
ratios for [Fe/H]<-2.5 imprinted by the nucleosynthe- 
sis in a few supernovae on the timescale ~ 10 7 yr and 
the other for [Fe/H]>-2 results from the mixing of the 
products from a whole site of the nucleosynthesis, taking 
place on the timescale longer than 10 9 yr. 
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